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I. INTRODUCTION 

As is well known, warm inflation -as opposed to the conventional "cool" inflation has 
the attractive feature of not necessitating a reheating phase at the end of the accelerated 
expansion thanks to the decay of the inflaton into radiation and particles during the slow-roll 
j^]. Thus, the temperature of the Universe does not drop dramatically and the Universe can 
smoothly proceed into the decelerated, radiation-dominated era essential for a successful 
big-bang nucleosynthesis [3|. This scenario has further advantages, namely: (i) the slow- roll 
condition <fi 2 <C V(<f>) can be satisfied for steeper potentials, (ii) the density perturbations 
originated by thermal fluctuations may be larger than those of quantum origin 4], and (iii) 
it may provide a very interesting mechanism for baryogenesis [5]. 

Warm inflation was criticized on the basis that the inflaton cannot decay during the 
slow-roll |6j. However, in the recent years, it has been demonstrated that the inflaton can 
indeed decay during the slow-roll phase -see ?J and references therein- whereby it now rests 
on solid theoretical grounds. 

Usually, for the sake of simplicity, when studying the dynamics of warm inflation the par- 
ticles created in the decay of the inflaton are treated as radiation thereby ignoring altogether 
the existence of particles with mass in the fluid thus generated. However, the very existence 
of these particles necessarily alters the dynamics as they modify the fluid pressure in two 
important ways: (i) its hydrodynamic, equilibrium, pressure is no longer p = p/3, with p 
the energy density of the matter-radiation fluid, but the slightly more general expression 
p = (7 — l)p where the adiabatic index, 7, is bounded by 1 < 7 < 2. (ii) It naturally arises 
a non-equilibrium, viscous, pressure n, via two different mechanisms: (a) the inter-particle 
interactions and (b) the decay of particles within the fluid {9]. 

A radiative fluid (e.g., a mixture of photons and electrons) is a well-known example of 
mechanism (a). Actually, it plays a significant role in the description of the matter-radiation 
decoupling in the standard cosmological model [3, 10]. Likewise, a hefty viscous pressure 
arises in mixtures of different particles species, or of identical species but with different 
energies -a case in point is the Maxwell-Boltzmann gas 111 ]. 

Concerning mechanism (ii), it is well known that the decay of particles within a fluid 
can be formally described by a bulk viscous pressure, n. This is so because the decay is 
an entropy-producing scalar phenomenon linked to the spontaneous widening of the phase 



2 



space and the bulk viscous pressure is also an scalar entropy-producing agent. In the case 
of warm inflation it has been proposed that the inflaton can excite a heavy field and trigger 
the decay of the latter into light fields Jjj]. 

Recently, a detailed analysis of the dynamics of warm inflation with viscous pressure 
showed that when II 7^ the inflationary region takes a larger portion of the phase space 
associated to the autonomous system of differential equations than otherwise [13J. It then 
follows that the viscous pressure facilitates inflation and lends support to the warm infla- 
tionary scenario. 

Our main target is to study the scalar and tensorial cosmological perturbations associated 
to this scenario and contrast them with the data gathered by the the three-year Wilkinson 



Microwave Anisotropy Probe (WMAP) 



14j . We shall follow the method used in a recent 



paper by two of us in which the cosmological perturbations generated by warm inflation 
driven by a tachyon field were considered 15( . The main differences between the latter work 
and the present one are: (i) here the inflaton is a scalar field, not a tachyon field; (ii) in Ref. 
[jjj no bulk dissipative pressure was considered. For the viscous pressure we shall assume 
the usual fluid dynamics expression II = —3(H js], where ( denotes the phenomenological 
coefficient of bulk viscosity and H the Hubble function. This coefficient is a positive-definite 
quantity (a restriction imposed by the second law of thermodynamics) and in general it is 
expected to depend on the energy density of the fluid. We shall resort to the WMAP data 
to restrict the aforesaid coefficient. 

The outline of the paper is as follows. Next Section presents our model. Sections III 
and IV deal with the scalar and tensor perturbations, respectively. Section V specifies to a 
chaotic potential and studies the limit of high dissipation. Finally, Section VI summarizes 
our findings. We choose units so that c = h = 8ttG = 1. 



II. WARM INFLATIONARY MODEL WITH VISCOUS PRESSURE 

We consider a spatially flat Friedmann- Robert son- Walker (FRW) universe filled with a 
self- interacting scalar field <p (the inflaton), of energy density, p^ = ^<fi 2 + V(<p), and an 
imperfect fluid, of energy density p and total pressure p + II, consisting into a mixture of 
matter and radiation of adiabatic index 7. 
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The corresponding Friedmann equation reads 

3# 2 =^ + V(<P)+p. 



(1) 



Because of the inflaton decays with rate T into the imperfect fluid the conservation equations 
for the inflaton and the fluid generalize to 

fa + m{ P4> + P4> ) = -r<p 2 4> + (3H + r)0 = -v >(f> , (2) 

and 

p + 3H(p + p + Tl) = p + 3H( 1P + U) = T0 2 , (3) 

respectively. In general, T depends on <p and as a consequence of the second law of ther- 
modynamics one has that T = T(0) > 0. As is customary, an upper dot indicates temporal 
derivative, and V<p = dV(<fi)/d<f). 

During the inflationary phase V(<f)) dominates over any other form of energy, therefore 
Friedmann's equation reduces to 

3H 2 = V{<j>) . (4) 

Likewise, imposing the slow-roll conditions, <p 2 V((p), and (3H + T)<j) [1], Eq. (T2J) 
becomes 

3F[l + r]0=-y^, (5) 

where the quantity r = T/(3H) quantifies the dissipation of the model. Throughout this 
paper we shall restrict our analysis to the high dissipation regime, i.e., r ^> 1. The reason for 
this limitation is the following. Outside this regime radiation and particles produced both 
by the decay of the inflaton and the decay of the heavy fields will be much dispersed by 
the inflationary expansion, whence they will have little chance to interact and give rise to a 
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non-negligible bulk viscosity. Likewise, because a much lower number of heavy fields will be 
excited the number of decays of heavy fields into lighter ones will diminish accordingly. (The 
weak dissipation reg ime (r < 1) has been considered by Berera and Fang (second reference 
of [2I) and Moss fie]]). Further, if r is not big, the fluid will be largely diluted and the mean 
free path of the particles will become comparable or even larger than the Hubble horizon. 
Hence, the regime will no longer be hydrodynamic but Knudsen's and the hydrodynamic 
expression n = —3(H we are using for the viscous pressure will become invalid. In such 
a situation, a consistent analysis should make use of Boltzmann's equation but this would 
complicate matters tremendously. This lies beyond the scope of this paper. 

Assuming that the decay of the inflaton is quasi-stable during inflation it follows that 
p <C 3if(7p + n) and p -C Tcj) 2 whence, from Eq. (jSj), we get 

rd) 2 - n 
p ~ 

7 



(6) 



for the energy density of the fluid. It becomes apparent that the viscous pressure (being 
necessarily negative) augments the energy density of the matter-radiation fluid. 
With the help of the dimensionless slow-roll parameter 



H 

IP 



2(1 



V 



1 2 



(7) 



a useful relation follows between the energy densities, namely, 



1 

P = - 

7 



2r 



n 



(8) 



.3(1 + r) r 

By definition, inflation lasts so long as a > 0. This amounts to the condition e < 1 
which with the help of last equation can be expressed as 



>3<I±lWn]. 



(9) 



Consequently, warm inflation with viscous pressure comes to a close when 



(1 + r) 

s^^bP + n]. (io) 



The number of e-folds at the end of inflation is given by 



N(<f>) 



/ _(l + r)#' ; 



(11) 



where the subscripts % and / stand for the beginning and end of inflation, respectively. 
For later purpose we introduce here the second slow-roll parameter, 77 = —777; 

, as a 

function of the potential and its two first derivatives, 



(1 + r) 



V, 



V 



I (Yd 

2\V 



(12) 



III. SCALAR PERTURBATIONS 



In terms of the longitudinal gauge, the perturbed FRW metric can be written as 



ds 2 = (1 + 2$) dt 2 - a{t) 2 (l - 2^) Sij dx i dx j 



(13) 



where the functions $ = $(£, x) and \l/ = x) denote the gauge- invariant variables 



of Bardeen 



17|. Introducing the Fourier components e , with k the wave number, the 



following set of equations, in the momentum space, follow from the perturbed Einstein field 
equations -to simplify the writing we omit the subscript k- 



$ + H$ = - 
2 



(jp + 11) a v 
k 



+ 



(14) 
(15) 



(6<f>) + [3H + T] (5(j)) + 



ru2 



,2 



5(f) = 40 $ - 



r + 2V, 



$, (16) 



(5p) + 3-fH5p + ka{ 1P + U)v + 3( 7 p + n)$ - 2 r,^0 - T0[2(50) + 0$] = 0, (17) 
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and 



v + AHv + 



8p 



+ 



T0 



(p + P) (P + P) 



5<p 



(18) 



where 



5 P = ( 7 - + <jn , 5n = n 



(19) 



and the quantity v arises upon splitting the velocity field as 5uj = ——n 1 v e lkx (J = 1, 2, 3) 



17|. 



Since the inflaton and the matter-radiation fluid interact with each other isocurvature 
(i.e., entropy) perturbations emerge alongside the adiabatic ones. This occurs because warm 



inflation can be understood as an inflationary model with two basics fields 18|, ll9| . In this 



context, dissipative effects themselves can produce a variety of spectral ranging from red 



to blue 
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20], thus producing the running blue to red spectral suggested by WMAP 
three-year data [14| . We will come back to this point below. 

When looking for non- decreasing adiabatic and isocurvature modes on large scales, 
k aH (which depend only weakly on time), it is permissible to neglect $ and those terms 
with two-times derivatives. Upon doing this and using the slow-roll conditions, the above 



equations simplify to 



2 H 



(jp + 11) a v 

k 



+ 



(20) 



[3H + r] (<tyj + 



<f>r + 2V, 



(21) 



Sp 



3-fH 



[r ^0 + r$], 



(22) 



and 



7 



AaH 



. + ,(7-i)fr + flT| + _rg =f , 



7p + n 



7p + n 



By combining of Eqs. (1221) and (123]) with (1201) . the latter becomes 



(23) 



where 



2HG((j)) 
G{<j>) = 1 



^ + (>- 1|+ii t)t^ 



8# 2 



2 7 p + 3n + 



7 



(24) 



(25) 



In the non-viscous limit, i.e., II — > when G(0) — > 1 (since 1 ^> p/H 2 ) and 7 = 4/3, 
Eq. (|24l reduces to the expression obtained for the conventional warm inflation 18]. 

With the help of Eq. (l2~41) and substituting by t as independent variable, Eq. (T2~Tj) can 
be solved. Upon using Eq. (jSJ) we get 



(26) 



and introducing the ancillary function 



V = 77- exp 



(3# + r) 



(27) 



we are led to the following differential equation for (p, 



_ 3 (r + 6H) 
tp 8 G{4>) (r + 3H) 2 _ 



r + 4H- I ( 7 -i) + n^ 



C / 3 1 H{3H + T) 



2* (28) 



Upon integrating it and resorting to Eq. ( 127ft we get 



where $s(4>) is given by 



3(0) 



(3# + r) 



C^exp[3(0)], 



3 (r + 6#) 



8 G{(j>) (r + 3#) 2 

(7-1 



T + AH 



3jH(3H + T)J V 



V , 



d t 



(29) 



(30) 
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and C is an integration constant. 

ThuMWe^pe—o. reads fl 



167T exp[— $5 



5 



V, 



(31) 



It is seen that in the absence of inflaton decay (T = 0, and therefore, II = 0), Eq. (IBTj) 
reduces to the typical expression for scalar perturbations, 5h ~ V5<f>/(H<j)) ~ H5<fi/<j). 
For high dissipation (r ^> 1), Eq. fl3T|) simplifies to 



167r\ 2 exp[— 2 3 



15 



r 2 



(32) 



where now 3(0) := 3(0) | r> i becomes 



3(0) = - 



1 



■r.«A + 



(7 



C / 97riJ 2 



3#r~"" ' 8G(0) 
During the slow- roll phase, the total density fluctuation, Spp 

perturbation, $, are related -at first approximation- by 



Qn(V))Ad<f>. (33) 
+ 5p, and the metric 



8p 



T 



V, 



-2[l+r]VG(0) 



1 + _L + ( [7 _ 1] + n ^)l^ 



(34) 



where use of Eq. (1241) has been made. Again, for r = IT = 0, we recover the usual relation 
Spr/pr — —2$, typical of cool inflation, in which pr = p<f> + p — V(<f>). 

In warm inflation, the fluctuations of the scalar field arise mainly from thermal interac- 



tion with the matter-radiation field. Therefore, following Taylor and Berera 
we can write 



2l|], for r > 1 



{5^ 



kp T r 
2tt 2 



(35) 



where T r stands for the temperature of the thermal bath and the wave number kp is given 
by kp = y/TH = H \/3r > H, and corresponds to the freeze-out scale at which dissipation 
damps out the thermally excited fluctuations. The freeze-out wave number kp is defined at 



the point where the inequality V aa < Eff, holds 



2l|. 



From Eqs. (151]) . (1521 and (1531) it follows that 



SI 



25tt : 



■ exp[— 23 



r r l/2 yzj2 



(36) 



where e ~ 77- 

2r 



Ik 

v 



denotes the dimensionless slow-roll parameter in the high dissipation 
phase. 

The scalar spectral index n s is given by 



n s - 1 



din 8 2 H 
dink 



(37) 



where the interval in wavenumber and the number of e-folds are related by 
dlnk(cf)) = —dN(<p). Upon using Eqs. fl36|) and fl3?j) . it can be written as 



7?,,. 



1 - 



e + 2rj + 



where 







_ r <<i>~ 








2r. 



(38) 



r] fa - 

r 



V, 



V 



1(Ya 

2 1 v 



(39) 



stands for the second slow- roll parameter, r/, when r ^> 1. 

One interesting feature of the three-year data gathered by the WMAP experiment is a 
significant running in the scalar spectral index dn s ldlnk = a s Ijj]. Dissipative effects can 
lead to a rich variety of spectral from red to blue [ji], [2^] . From Eq. (1381) it is seen that in 
our model the running of the scalar spectral index is given by 



a. 



'2e ... 
r 



£ 

r 



23 s- — 

10 2r 



+ 



43 



(ln(r)), 



• (40) 



e r 

In models with only scalar fluctuations, the marginalized value of the derivative of the 
spectral index can be approximated by dn s /dlnk = a s ~ —0.05 for WMAP-3 only [l4j . 
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IV. TENSOR PERTURBATIONS 



As argued by Bhattacharya et al. 22j, the generation of tensor perturbations during 
inflation gives rise to stimulated emission in the thermal background of gravitational waves. 
As a consequence, an extra temperature dependent factor, coth(/c/2T), enters the spectrum, 



A 2 oc k ns . Thus it now reads, 



A 





2 


" k ' 


V 


" k ' 




| coth 


~ — — coth 




_2T_ 


O 7T 2 


_2T_ 



(41) 



the spectral index being 



d 



d Ink 



In 



A 2 

A 9 



coth[A;/2T] 



-2e 



(42) 



where we have used Eq.(J7J). 

A quantity of prime interest is the tensor-scalar ratio, defined as R(ko) = (^j^ 
where Pr = 255^/4 and k is known as the pivot point. Its expression in the high 
dissipation limit, r ^> 1, follows from using Eqs. (1361) and (14ip . 



k=ko 



R(k ) 



A 2 



P', 



n 



k=kg 



r r l/2 y5/2 



Combining the WM AP three-year data 



exp[2 $S((j))] coth 



k 

2T 



(43) 



k=ko 



141 ] with the SDSS large scale structure surveys 



231 ] . yields the upper bound i?(/;;o=0.002 Mpc _1 )< 0.28 at 95% confidence level, where 
ko = 0.002 Mpc -1 corresponds to r k ~ 30, with the distance to the decoupling surface 
being tq = 14.4 Mpc. The SDSS is a measure of the galaxy distributions at red-shifts 
a ~ 0.1 and probes k in the range 0.016 h Mpc _1 < k < 0.011 /iMpc" 1 . The recent WMAP 
three-year data imply for the scalar curvature spectrum the values Pn(ko) = 25 5^(&;o)/4 ~ 
2.3 x 10~ 9 and the tensor-scalar ratio R(ko) = 0.095. We shall make use of these values to 
set constraints on the parameters of our model. 
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V. THE HIGH DISSIPATION LIMIT WITH CHAOTIC POTENTIAL 



In this section we assume for the inflaton the chaotic potential 



V(0) = -m 2 2 , 



(44) 



where m > is a free parameter, and (as mentioned above) we restrict ourselves to study 
the high dissipation regime (r ^> 1). 



A. The case r = T , C = Co 

When the inflaton decay rate and the bulk viscosity coefficient are given by T = T 
constant > and ( = (q= constant, respectively, the first slow roll parameter reduces to 



v^6 m 1 



Likewise, the Hubble factor is given by 



and the parameter r becomes 



H{<j>) 



m 4> 



o in d) 



(45) 



(46) 



(47) 



Because in this scenario 



3rH 



(48) 



the inflaton depends on time as 

(pit) = 4>i exp 



m 2 1 




r ii 

m 






1 t 


[ r . 




r 



(49) 



where <pi = <fr(t = = 0). This entirely agrees with Taylor and Berera findings 21 1 
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As for the energy density of the matter-radiation fluid we have the expression 



m 4> 

7 



\/E m 2 

3r n 



+ 



Ho 



(50) 



or in terms of ps, 



P 



V3 

7 



2 m 2 

~3r\T 



+ 3Co 



1/2 



(51) 



By means of Eq. ( TIT]) , the number of e-folds at the end of warm inflation is found to be 



total 



<t>f y 



3 m 



(52) 



where, because of > V/, the initial and final values of the scalar field satisfy 0, > <pf. 

It is readily seen that e/ ~ 1 (i.e., a(i = £/) ~ 0), and that the scalar field at the end of 
inflation reads 



m 



r 



(53) 



Rewriting the total number of e-folds in terms of fa and 0/ and using Eq. (l53|) . we get 



-[N total + 2] <f> f 



a/6 m r 

A 1 o 



(54) 



To comfortably solve the problems of the standard big-bang cosmology we must have N to tai 
60 e-folds. This implies fa ~ 76m/r . 
At the beginning of inflation 



r(<j ) = fa) = r i = ±- (^] > I. 

93 V m 



(55) 



which implies that To ^> v93 m. 
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From Eq. (!36|) . the scalar power spectrum results to be 



Pn(k ) » ^- 2 [^T V{^y' 2 + 2V3m 2 (l - 2 7 ) + 3^3^(2 - 3 7 )] 3/2 



pl/2 

1 Q J r 



3V4 m 2 y (0 Q )3/4 



(56) 



Likewise, Eg. (143 j) provides us with the tensor-scalar ratio 



R(ko) ~ -[8 7 r o ^(0o) 1/2 + 2V3m 2 (l - 2 7 ) + 3^0^(2 - 3 7 )]" 



-3/2 



X 



3 1 / 4 mV(0 o ) 7/4 



-pl/2rp 
1 Q J r 



coth 



2T 



(57) 



where V((po) and 0o stand for the potential and the scalar field, respectively, when the 
perturbation, of scale ko = 0.002Mpc _1 , was leaving the horizon. 

By resorting to the WMAP three-year data, Pn{.h) ^ 2.3 x 1(T 9 and R(k ) = 0.095, 
and choosing the parameters 7 = 1.5, m = 10 _6 mp, T ~ T r ~ 0.24 x 10 16 GeV and 
k = 0.002 Mpc~\ it follows from Eqs. ([56]) and ([57} that V((j) Q ) ~ 1.5 x lO^mf, and 
( ~3x 10 _6 mp. When the scale ko was leaving the horizon the inflaton decay rate r is 
seen to be of the order of 10~ 3 mp. Thus Eq. ( 140]) tells us that one must augment Co by two 
orders of magnitude to have a running spectral index a s close to the observed value 14]. 



B. The case T = r(</>), C = C(p) 

Here we assume, T = T(0) = a V(4>) = am 2 2 /2, and ( = ((p) = Ci p, where a and (j 
are positive-definite constants. Obviously, II = —3(ip H. 

Using the chaotic potential, Eq. (|44p . we find that the slow roll parameter, the Hubble 
factor and r can be written as 

_ 2\/6 „, md) , mad) _ 

e=—^— H(#) = -%, and r = > 1 , 58 

respectively. By combining the middle equation of ([58]) with Eq. (j4~8]) we find that 



2 (t) = 2 
14 



(59) 



The energy density of the matter-radiation fluid reads 



P 



4v/3 



m 



30 [v^ 7 -3Cim< 



(60) 



and, in terms of p^, 



12 m 2 



3p 



1/2 



1 



T-3& 



1/2 



(61) 



Using Eq. pip , the number of e- folds at the end of warm inflation results to be, 



total 



4>f v 



am 



r 36 

For 6f ~ 1 it is seen that the scalar field at the end of inflation reduces to 



(62) 



2VQ 



a m 



1/3 



(63) 



Rewriting the total number of e- folds in terms of <ftf and (pi, and using last equation, we 
obtain 



[3AU^ + 1] 1/3 <p f = [3N total + l] 



1/3 



2^ 



1/3 



(64) 



To get N tota i « 60 e-folds we must have ^ « [362^/(0 m)] 1 / 3 . 
At the beginning of inflation the r parameter becomes 



r((f) = (pi) = Ti 



181 



nl/3 



a 2 / 3 m 2 / 3 



(65) 



resulting in the requirement that a ^> y 181m i m the high dissipation regime. 
From Eq. fl36|) . the scalar power spectrum is seen to obey 



Pn(ko) 



2tt 2 



3?n 



a 



1/2 T 



(66) 



with 



% = o? V(0 O ) 2 + ^*» a 6)/(«Ci-*r) + ^ + |y _ ^ _ |y j 
being d 7^ 87/9, and where 

ft = -87^0) + 6 7 2 - 12^7(1^0)- 3 7 + 3V3{iVWo) + 
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fa = 2v^^)(8 7 - 9d 2 ) + 3V3Ci(4 7 - 1), 
fa = (8 7 - 9d 2 ) yj 192 7 3 + 216 7 2 Ci 2 " 96 7 2 - IO87C1 2 + ^Cf , 

and 

(3 A = I ] [64m 2 7 2 + 27m 2 7 C 2 - 64m 2 7 + 36m 2 Ci + 72(i 7 2 - 18& 7]. 
V 3p3 J 

From Eq. (l43p the tensor-scalar ratio can be obtained as 



R(k ) » I 



3 1 / 4 m 2 F(0 O ) 5/4 



coth ( A ) . (67 ) 



aV2 K 7; 

Resorting again to the WMAP three-year data (P n (h) ^ 2.3 x 10~ 9 , R(k ) = 0.095) and 
choosing the parameters 7 = 1.5, a ~ 10 7 mp 3 , V(4>o) ~ 1.5xl0 _11 mp, T ~ T r ~ 0.24xl0 14 
GeV and k = 0.002Mpc _1 , we find from Eqs. fl66|) and fl67|) that m ~ 10 3 mp, and 
Ci ~ 10 _8 mp 1 . We note that the dissipation coefficient when the scale k$ was leaving the 
horizon is of the order of r(0 o ) = a V(<fro) ~ 10~ 4 mp. It follows from Eq. (1401) that one must 
decrease Ci by three orders of magnitude to have a running spectral index a s close to the 
WMAP observed value. 

VI. CONCLUDING REMARKS 

In this paper we considered a warm inflationary scenario in which a viscous pressure is 
present in the matter-radiation fluid. This pressure arises on very general grounds, either as 
a hydrodynamical effect 8] or as a consequence of the decay of massive fields -previously ex- 
cited by the inflaton- into light fields, or both. We investigated the corresponding scalar and 
tensor perturbations. The contributions of the adiabatic and entropy modes were obtained 
explicitly. Specifically, a general relation for the density perturbations is given in Eq. (!3TT) . 
The tensor perturbations are generated via stimulated emission into the existing thermal 
background, Eq. (14T1) . and the tensor-scalar ratio -as well as the dissipation parameter- is 
modified by a temperature dependent factor. 

The effect of the viscous pressure reveals itself at the e-folding level. Indeed, in the 
first case, i.e., when T = r = constant and ( = (q, as Eq. (l52|) shows, the total number 
of e-folds depends on the difference of the inflaton field evaluated at the beginning and 
at the end of inflation, i.e., <bj — <bf, where the initial value of the inflaton field is 6"} = 



16 



^P^(3Hi + r ) 2 (7 pi + ILJ. Because IT < one has that (pi is lower than when the viscous 



pressure is not considered, thereby the viscous pressure shortens the inflationary phase. 
Also, in this first case, we constrained the parameters of our model with the help of the 



WMAP three-year data 14j . Thus, we have found that £q should be the order of 10~ 6 nip, 
and the potential, V((f>o), of the order of 10 -11 mf, when the perturbation exits the horizon at 
the scale of ko = 0.002 Mpc^ 1 . Likewise, we have found that the contribution of the viscous 
pressure to the scalar power spectrum, i.e., x = — j is the order of two percent. 

Similar results follow in the second case, i.e., when T = T(0) and C = C(p)- 111 this 
instance, the ratio \ is about four percent. While the deviation in both cases is small (but 
not negligible), we think it is still big enough to be detected by future measurements of the 
large scale structure of the Universe. 

We have implicitly assumed that 7 is a constant which is not really the case since matter 
and radiation redshift at different rates, and that the bulk viscosity pressure is given by the 
conventional stationary formula U = —3(H while it is well known that it does not respect 
relativistic causality and should be generalized to the Israel-Stewart expression rll + IT = 



—3(H, where r is the relaxation time of the viscous process 24]. However, it is to be 
expected that these generalizations, when properly incorporated, will not substantially alter 
our findings. 

We have not studied the effect of the viscous pressure on the bispectrum of density 
perturbations -though there are good reasons to do it- as it lies a bit outside the main target 
of our paper. While cool inflation typically predicts a nearly vanishing bispectrum, and 
hence a small (just a few per cent) deviation from Gaussianity in density fluctuations -see 
e.g. (25]-, warm inflation clearly predicts a non- vanishing bispectrum. The latter effect arises 
from the non-linear coupling between the the fluctuations of the inflaton and those of the 
radiation. This can produce a moderate non- Gaussianity [26( or even a stronger one -likely 
to be detected by the Planck satellite (23]- if the aforesaid nonlinear coupling is extended 



to subhorizon scales 



28( | . Because II implies an additional coupling between the radiation 



and density fluctuations it is to be expected that non- Gaussianity will be further enhanced. 
Perhaps, this could serve to observationally constrain II by future experiments. We plan to 
consider this question in a subsequent paper. 

In summary, our model presents two interesting features: (i) Related to the fact that the 
dissipative effects plays a crucial role in producing the entropy mode, they can themselves 
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produce a rich variety of spectral ranging from red to blue. The possibility of a spectrum 
which does run so is particularly interesting because it is not commonly seen in inflationary 
models which typically predict red spectral, (ii) The viscous pressure may tell us about 
how the matter-radiation component behaves during warm inflation. Specifically, it will 
be very interesting to know how the viscosity contributes to the large scale structure of 
the Universe. In this respect, we anticipate that the Planck mission 27| will significantly 
enhance our understanding of the large scale structure by providing us with high quality 
measurements of the fundamental power spectrum over an larger wavelength range than the 
WMAP experiment. 



Acknowledgments 

SdC and RH wish to thank the "Universidad Autonoma de Barcelona" for hospitality 
where part of this work was done. SdC was supported by Comision Nacional de Ciencias y 
Tecnologfa through FONDECYT grants N° 1030469, and "FONDECYT-Concurso incentivo 
a la cooperacion internacional" N° 7060005. Also, from UCV-DGIP N° 123.764. RH was 
supported by the "Programa Bicentenario de Ciencia y Tecnologfa" through the Grant 
"Insertion de Investigadores Postdoctorales en la Academia" N° PSD/06. DP acknowledge 
"FONDECYT-Concurso incentivo a la cooperacion internacional" N°7060005, and is grateful 
to the "Institute de Fisica" de la PUCV for warm hospitality, also DP research's was partially 
supported by the "Ministerio Espanol de Education y Ciencia" under Grant FIS2006- 12296- 
C02-01. 



[1] E.W. Kolb and M.S. Tuner, The Early Universe (Addison- Wesley, Reading, MA, 1990); D. 
Lyth, Cosmological Inflation and Large Scale Structure (Cambridge University Press, Cam- 
bridge, UK, 2000); A. Linde, jhep-th/0503203[ 

[2] A. Berera, Phys. Rev. Lett. 75, 3218 (1995); A. Berera and L.Z. Fang, Phys. Rev. Lett. 74, 
1912 (1995); A. Berera, Phys. Rev. D 55, 3346 (1997). 

[3] P.J.E. Peebles, Principles of Physical Cosmology (Pricenton University Press, Princeton, 
1993). 



18 



[4] A. Berera, Nucl. Phys. B585, 666 (2000); A.N. Taylor and A. Berera, Phys. Rev. D 62, 083517 
(2000); S. Gupta, A. Berera, A.F. Heavens, and S. Matarrese, Phys. Rev D 66, 043510 (2002). 
[5] R. H. Brandenberger and M. Yamaguchi, Phys. Rev. D 68, 023505 (2003). 
[6] J. Yokoyama and A. Linde, Phys. Rev. D 60, 083509 (1999). 

[7] A. Berera and O. Ramos, Phys. Rev. D 71, 023513 (2005); L.M. Hall and I.G. Moss, Phys. 

Rev. D 71, 023514 (2005); M. Bastero-Gil and A. Berera, Phys. Rev. D 71, 063515 (2005). 
[8] L. Landau and E.M. Lifshitz, Mecanique des Fluides (MIR, Moscow, 1971); K. Huang, Sta- 
tistical Mechanics (J. Wiley, New York, 1987). 
[9] Ya. B. Zeldovich, JEPT Lett. 12, 307 (1970); J.D. Barrow, Nucl. Phys. B310, 743 (1988); 
W. Zimdahl and D. Pavon, Phys. Lett. A 176, 57 (1993); W. Zimdahl, Mon. Not. R. Astron. 
Soc. 280, 1239 (1996). 
[10] S. Weinberg, Gravitation and Cosmology (J. Wiley, New York, 1972), pp. 51-58. 
[11] S. Harris, An Introduction to the Theory of Boltzmann Equation (Holt, Reinhart and Winston, 
New York, 1971); C. Cercignani, Theory and Applications of the Boltzmann Equation (Scotish 
Academic Press, Edinburgh, 1975). 
[12] A. Berera and R.O. Ramos, Phys. Lett. B 567, 294 (2003); Phys. Rev. D 71, 023513 (2005). 
[13] J.P Mimoso, A. Nunes, and D. Pavon, Phys. Rev. D 73, 023502 (2006). 
[14] D.N. Spergel, et al, |astro-ph/0603449 

[15] R. Herrera, S. del Campo and C. Campuzano, JCAP 10(2006)009. 
[16] I. G. Moss, Phys. Lett. B 154, 120 (1985). 
[17] J. Bardeen, Phys. Rev. D 22, 1882 (1980). 

[18] A. Starobinsky and J. Yokoyama, in Proceedings of the Fourth Workshop on General Relativity 
and Gravitation, edited by N. Nakao et al. (Kyoto University, Kyoto, 1995) pp. 381; A. 
Starobinsky and S. Tsujikawa, Nucl. Phys. B 610, 383 (2001); H. Oliveira and S. Joras, Phys. 
Rev. D 64, 063513 (2001). 

[19] H. Oliveira, Phys. Lett. B 526, 1 (2002). 

[20] L.M.H. Hall, I.G. Moss, and A. Berera, Phys. Rev. D 69, 083525 (2004); I.G. Moss, Phys. 

Lett. B 154, 120 (1985); 
[21] A.N. Taylor and A. Berera, Phys. Rev. D 55, 3346 (1997). 

[22] K. Bhattacharya, S. Mohanty, and A. Nautiyal, Phys. Rev. Lett. 97, 251301 (2006). 
[23] M. Tegmark, et al, Phys. Rev. D 69, 103501 (2004). 



19 



[24] W. Israel, Ann. Phys. (NY) 100, 310 (1976); W. Israel and J.M. Stewart, Ann. Phys. (NY) 
118, 341 (1979); D. Pavon, D. Jou, and J. Casas- Vazquez, Ann. Institute H. Poincare A 36, 
79 (1982). 

[25] A. Gangui, F. Lucchin, S. Matarrese, and S. Mollerach, Astrophys. J. 430, 447 (1994). 

[26] S. Gupta, A. Berera, A.F. Heavens, and S. Matarrese, Phys. Rev. D 66, 043510 (2002); S. 

Gupta, Phys. Rev. D 73, 083514 (2006). 
[27] |http://www.rssd.esa.int/SA/PLANCK/docs/JASR2004.pdf 
[28] I.G. Moss, and Ch. Xiong, |astro-ph/070~1302j 



20 



